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Abstract We describe Odorrana liboensis sp. nov., a new 
species from the Maolan National Nature Reserve, Libo 
County, Guizhou Province, China. Phylogenetic analyses 
based on DNA sequences of the mitochondrial 12S rRNA, 
16S rRNA, and ND2 genes supported the new species as 
an independent lineage, closely related to O. lipuensis. The 
uncorrected genetic distances between the 12S rRNA and 
16S rRNA in the new species and its closest congener, O. 
lipuensis, were 6.06% and 5.19%, respectively. The new species 
is distinguished from its congeners by a combination of 
the following morphological characters: (1) having medium 
body size, with the snout-vent length (SVL) of adult 
females approximately 1.2 times as long as males at 56.9 + 
10 (55.8-58.2 mm, n = 9) in females and 48.7 + 12 (471-499 
mm, n = 5) in males; (2) head length greater than width in 
males and females; (3) tympanum distinctly visible, greater 
than one-half the diameter of the eye; (4) eyes big and 
prominent, width of upper eyelid (UEW) approximately 
3/4 of eye diameter (ED); (5) dorsolateral folds absent; (6) 
two metacarpal tubercles; (7) relative finger lengths: II < I < 
IV < III; (8) subarticular tubercles on fingers prominent: 1, 1, 
2, 2; (9) one metatarsal tubercle; (10) tibiotarsal articulation 
reaching to between the eye and nostril when the leg is 
stretched forward; (11) toes with entire webbing to disks; 
(12) subarticular tubercles on toes prominent: 1, 1, 2, 3, 2; (13) 
dorsal surfaces of limbs with distinct brownish-black bands; 
(14) smooth, grass-green dorsum with irregular brown 
mottling; (15) venter smooth, lacking black spots; and (16) 
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lacking pectoral spinules, lacking vocal sacs, and light white 
nuptial pad present on finger I in males. The new species is 
currently only known from the type locality. 


Keywords Taxonomy, morphology, Odorrana, Odorrana 
liboensis sp. nov., karst cave, Guizhou 


1. Introduction 


The genus Odorrana was first recognized by Fei et al. (1990) 
with the type species Odorrana margaretae (Liu, 1950). Despite 
previous controversy over classification (Dubois, 1992; Matsui, 
1994), molecular phylogenetic analyses suggested that Odorrana 
forms a well-supported monophyletic group (Matsui et al., 
2005; Ngo et al., 2006; Cai et al, 2007; Che et al., 2007; Stuart, 
2008; Chen et al, 2013; Li et al, 2018a). According to Frost (2021), 
Odorrana contains at least 61 recognized species and is widely 
distributed in montane streams and rivers in the subtropical 
and tropical regions of East and Southeast Asia (Fei et a, 2012; 
AmphibiaChina, 2021; Frost, 2021). The range of species in 
this genus covers the Ryukyu Archipelago, southern China, 
northeastern India, and the Thai-Malay Peninsula, and extends 
southwards to Sumatra and Borneo (Frost, 2021). All species in 
the genus are associated with mountain streams except for O. 
wuchuanensis (Xu, 1983) and O. lipuensis Mo, Chen, Wu, Zhang & 
Zhou, 2015, which occur in dark caves. 

Among the Odorrana species, 39 occur in China and 24 of 
these are endemic to China (Fei et al, 2012; AmphibiaChina, 
2021; Frost, 2021; Chen et al., 2020; Shen et al., 2020; Liu et 
al., 2021). Fei et al. (1990) established Odorrana, and later Ye 
and Fei (2001) suggested four species groups (O. andersonii, 
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O. kuangwuensis, O. schmacheri, and O. livida) within Odorrana 
based on a morphological phylogenetic study. Fei et al. (2005) 
established two subgenera (Odorrana and Bamburana) within 
Odorrana and recognized Odorrana versabilis as the type species 
of Bamburana. Molecular phylogenetic studies support the 
division of species groups within Odorrana but not the two 
subgenera (Che et al, 2007). Subsequently, Fei et al. (2009) divided 
the Chinese Odorrana species into six species groups (O. tormota, 
O. taiwaniana, O. graminea, O. margaretae, O. schmackeri, and O. 
andersonii; O. taiwaniana is a synonym for O. swinhoana) based on 
morphological characters. These divisions have been accepted by 
other researchers and applied to the Odorrana species distributed 
outside China (Pham et al., 2016; Li et al, 2018a). However, recent 
studies have rejected the monophyly of the O. margaretae, O. 
schmackeri, and O. andersonii species groups (Chen et al., 2013). 
The phylogenetic relationships between these species groups, 
the species included in the species groups, and the species group 
affiliations of new species published in recent years are unclear. 

The montane river and stream habitats of most Odorrana 
may promote intraspecific divergence through geographical 
isolation, especially given the wide geographical range of the 
complex (Chen et al., 2013; Li et al, 2015). Indeed, many cryptic 
species of Odorrana have been discovered. For example, 25 new 
species of Odorrana have been described since 2005 (See species 
list of the genus Odorrana in Frost (2021) and AmphibiaChina 
(2021)). This indicates that further discoveries or potential 
taxonomic changes within Odorrana are likely. 

From 2016 to 2018, 14 specimens of an unknown species of 
Odorrana were collected inside a completely dark karst cave in 
the Maolan National Nature Reserve, Libo County, Guizhou 
Province, China. Morphologically, these specimens most closely 
resemble O. lipuensis and O. kweichowensis Li, Xu, Lv, Jiang, 
Wei & Wang, 2018 (Mo et a, 2015; Li et al, 2018a), but differ 
from O. lipuensis, O. kweichowensis, and all other Odorrana from 
China and adjoining countries. They inhabit a dark karst cave 
environment similar to that of O. wuchuanensis and O. lipuensis 
but are morphologically very different from O. wuchuanensis 
and O. lipuensis. To distinguish these specimens, we conducted 
phylogenetic analyses based on mitochondrial DNA and 
morphological comparisons. All of the analyses consistently 
indicated that the specimens from Maolan National Nature 
Reserve are a new taxon. We describe this taxon here as a new 


species. 
2. Materials and Methods 


2.1. Sampling A total of 72 specimens were collected in this 
study. Fourteen specimens were of the undescribed species from 
Maolan National Nature Reserve in Libo County, Guizhou 
Province; 35 were O. kweichowensis from Lengshuihe Nature 


| Vol. 12 


Reserve in Jinsha County, Guizhou Province; four were O. 
yizhangensis Fei, Ye & Jiang, 2007 from Kuankoushui National 
Nature Reserve in Suiyang County, Guizhou Province, and 
Yuntai Mountain, Shibing County; one was O. huanggangensis 
Chen, Zhou & Zheng, 2010 from Yueliangshan Nature 
Reserve in Conjiang County, Guizhou Province; one was O. 
kweichowensis from Maolan National Nature Reserve in Libo 
County, Guizhou Province; and 17 were O. wuchuanensis from 
Maolan National Nature Reserve, Libo County, Guizhou 
Province (Figure 1). All of the specimens were fixed in 10% 
buffered formalin and later transferred to 75% ethanol for 
preservation. The muscles used for molecular analysis were 
preserved in 95% alcohol at —20 °C. All of the specimens were 
deposited in Guizhou Normal University (GZNU), Guiyang 
City, Guizhou Province, China. 


2.2. DNA Extraction, PCR amplification, and 
sequencing Genomic DNA was extracted from muscular 
tissue using a DNA extraction kit from Tiangen Biotech 
Co., Ltd. (Beijing). All samples were sequenced for three 
mitochondrial genes: partial 12S ribosomal RNA gene (12S 
rRNA), 16S ribosomal RNA gene (16S rRNA), and NADH 
dehydrogenase subunit 2 (ND2). The primers used for 12S 
rRNA were 12SF (5-GGTTTGRTCCTRGCCTTAC-3) 
and 12SR (5-CCATGTTACGACTTGCCTCT-3) following 
Chen et al. (2013), the primers used for 16S rRNA were 
16SF (5-ACGAGCCTAGTGATAGCTGGTT-3) and 
16SR (5-CGGTCTGAACTCAGATCACGT-3) following 
Chen et al. (2013), and the primers used for ND2 were Gln- 
LND2 (5-CCCTTTGCACTTCCTTTATGC-3) and Ala- 
HND2 (5-GGCCTGAGTTGCATTCATG-2) following 
Li et al. (2015). PCR amplifications were performed in a 25 pL 
reaction volume with the following cycling conditions: an initial 
denaturing step at 95 °C for five min; 36 cycles of denaturing 
at 95 °C for 40 s, annealing at 52 °C (for 12S rRNA)/49 °C (for 
16S rRNA)/60 °C (for ND2) for 40 s and extending at 72 °C for 
1 min, and a final extending step of 72 °C for 10 min. The PCR 
products were sequenced on an ABI Prism 3730 automated DNA 
sequencer in Chengdu TSING KE Biological Technology Co. Ltd. 
(Chengdu, China). All of the sequences have been deposited in 
GenBank (Table 1). Some homologous DNA sequences of 
voucher specimens of related species were downloaded 
from GenBank and incorporated into the phylogenetic 
analyses. 


2.3. Phylogenetic analyses We used a total of 187 
sequences (including 72 12S rRNA, 79 16S rRNA, and 36 ND2 
gene sequences) for molecular analyses. Three mitochondrial 
genes from 11 muscle tissues were sequenced and 154 sequences 
downloaded from GenBank from 48 species of the genus 
Odorrana were used. These included the undescribed species 
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from China that was named in this study (Figure 1). A total of 
26 seguences were downloaded from GenBank as out-groups 
(Chen et al., 2013). Detailed information on these materials is 
shown in Table 1. 

All of the sequences were assembled and aligned using the 
MUSCLE (Edgar, 2004) module in MEGA 7.0 (Kumar et al, 2016) 
with default settings. Alignments were checked by eye and revised 
manually if necessary. Trimming with the gaps partially deleted 
was performed in GBLOCKS 091b (Castresana, 2000). Phylogenetic 
trees were constructed with both maximum likelihood (ML) and 
Bayesian inference (BI). The ML was conducted in IQ-TREE 
(Nguyen et al, 2015) with 2000 ultrafast bootstrapping (Hoang et al, 
2018) and was performed until a correlation coefficient of at least 
099 was reached. The BI was performed in MrBayes 32.1 (Ronquist 
et al, 2012), and the best-fit model was obtained by the Bayesian 
inference criteria (BIC) computed with PartitionFinder 2 (Lanfear 
et al, 2016). For this analysis, DS rRNA, 16S rRNA, and ND2 genes 
were defined. 


Provincial boundary 
National boundaries 
Rivers 


Elevation (m) 
7411 
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The analysis suggested that the best partition scheme was 12S 
rRNA/I6S rRNA/ND2 genes. We selected GTR+I+G as the best 
model for 12S rRNA and 16S rRNA and the TIM+I+G model 
as the best model for the ND2 gene. Two independent runs 
were conducted in BI analysis, each of which was performed for 
2,000,000 generations and sampled every 1000 generations. The 
first 25% of the samples were discarded as burn-in, resulting in 
a potential scale reduction factor (PSRF) of «001. Nodes in the 
trees were considered well supported when Bayesian posterior 
probabilities (BPP) were 2095 and the ML ultrafast bootstrap value 
(UFB) was 29596. Uncorrected p-distances (1000 replicates) based on 
DSrRNA and 16S rRNA were calculated in MEGA 70 (Kumar et 
al, 2016). 


2.4. Morphological analysis Morphometric data were 
taken from 74 well-preserved adult specimens (voucher 
information in Table 2 and Table SI). Measurements were 


recorded to the nearest 0.1 mm with digital calipers by Tao 


e Odorrana liboensis sp. noy. 


@ Odorrana lipuensis ! 
@ Odorrana wuchuanensis 
PS Odorrana kweichowensis 
I Odorrana huakggangensis 
P Odorrana yizhangensis 


Figure 1 Sampling collection localities and distribution of the Odorrana liboensis sp. nov., O. lipuensis, O. wuchuanensis, O. kweichowensis, O. huanggangensis, 
and O. yizhangensis in southwest China. 1. Maolan National Nature Reserve, Libo County, Guizhou, China. 2. Lipu County, Guangxi, China. 3. Jiarong 
Town, Libo County, Guizhou, China. 4. Maolan National Nature Reserve, Libo County, Guizhou, China. 5. Lengshuihe Nature Reserve, Jinsha County, 
Guizhou, China. 6. Yueliangshan Nature Reserve, Congjiang County, Guizhou, China. 7. Yuntai Mountains, Shibing County, Guizhou, China. 8. Kuan- 
kuoshui National Nature Reserve, Suiyang County, Guizhou, China. The base maps are from Standard Map Service website (http://bzdt.chmnr.gov.cn/ 


index html). 
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Table 1 Localities, voucher information, and GenBank numbers for all samples used in this study. 


ID Species Locality Voucher 12S 16S ND2 


1 Odorrana nasuta 

2 Odorrana versabilis 

3 Odorrana exiliversabilis 
4 Odorrana nasica 

5 Odorrana yentuensis 
6 Odorrana tormota 

7 Odorrana amamiensis 
8 Odorrana narina 

9 Odorrana supranarina 
10 Odorrana swinhoana 


11 Odorrana utsunomiyaorum 


12 Odorrana aureola 

13 Odorrana livida 

14 Odorrana leporipes 

15 Odorrana chloronota 

16 Odorrana graminea 

17 Odorrana hosii 

18 Odorrana banaorum 

19 Odorrana morafkai 

20 Odorrana kweichowensis 
21 Odorrana kweichowensis 
22 Odorrana schmackeri 
23 Odorrana schmackeri 
24 Odorrana bacboensis 

25 Odorrana tiannanensis 
26 Odorrana fengkaiensis 
27 Odorrana hainanensis 
28 Odorrana nanjiangensis 
29 Odorrana hejiangensis 
30 Odorrana huanggangensis 
31 Odorrana huanggangensis 
32 Odorrana huanggangensis 
33 Odorrana tianmuii 

34 Odorrana grahami 

35 Odorrana junlianensis 
36 Odorrana andersonii 

37 Odorrana kuangwuensis 
38 Odorrana margaretae 
39 Odorrana jingdongensis 
40 Odorrana daorum 

41 Odorrana hmongorum 
42 Odorrana wuchuanensis 
43 Odorrana wuchuanensis 
44 Odorrana wuchuanensis 
45 Odorrana mutschmanni 
46 Odorrana yizhangensis 
47 Odorrana yizhangensis 
48 Odorrana yizhangensis 
49 Odorrana yizhangensis 
50 Odorrana yizhangensis 
51 Odorrana yizhangensis 
52 Odorrana lungshengensis 
53 Odorrana lungshengensis 
54 Odorrana anlungensis 
55 Odorrana chapaensis 

56 Odorrana geminata 

57 Odorrana ishikawae 

58 Odorrana absita 

59 Odorrana lipuensis 

60 Odorrana lipuensis 

61 Odorrana lipuensis 

62 Odorrana lipuensis 


63 Odorrana liboensis sp. nov. 
64 Odorrana liboensis sp. nov. 
65 Odorrana liboensis sp. nov. 


66 Amolops loloensis 

67 Amolops mantzorum 
68 Amolops granulosus 

69 Amolops sp. 

70 Amolops ricketti 

71 Amolops wuyiensis 

72 Amolops hongkongensis 
73 Hylarana guentheri 
74 Hylarana spinulosa 

75 Glandirana tientaiensis 
76 Pelophylax nigromaculata 
77 Babina adenopleura 

78 Babina daunchina 


Wuzhishan City, Hainan, China 

Leigongshan Nature Reserve, Leishan County, Guizhou, China 
Wuyishan City, Fujian, China 

HaTinh, Vietnam 

Guangxi, China 

Huangshan City, Anhui, China 

Tokunoshima, Ryukyu 

Okinawa Island, Japan 

Iriomotejima, Ryukyu 

Nantou County, Taiwan, China 

Iriomotejima, Ryukyu 

Phu Rua District, Loei Prov., Thailand 

Prachuap Kirikhan Prov, Thailand 

Shaoguan City, Guangdong, China 

Ha Giang, Vietnam 

Wuzhishan City, Hainan, China 

Kuala Lumpur, Malaysia 

Tram Lap, Vietnam 

TramLap, Vietnam 

Maolan National Nature Reserve, Libo County, Guizhou, China 
Lengshuihe Nature Reserve, Jinsha County, Guizhou, China 
Yichang City, Hubei, China 

Badagongshan Nature Reserve, Sangzhi County, Hunan, China 
Khe Moi, Nghe An, Vietnam 

Hekou County, Yunnan, China 

Heishiding Nature Reserve, Fengkai County, Guangdong, China 
Wuzhishan City, Hainan, China 

Yichang County, Hubei, China 

Hejiang County, Sichuan, China 

Yueliangshan Nature Reserve, Congjiang County, Guizhou, China 
Fanjingshan Nature Reserve, Jiangkou County, Guizhou, China 
Wuyishan Nature Reserve, Fujian, China 

Linn area , Zhejiang , China 

Kunming City, Yunnan, China 

Junlian, Sichuan, China 

Longchuan County, Yunnan, China 

Nanjiang County, Sichuan, China 

Emei City, Sichuan, China 

Jingdong County, Yunnan, China 

Sa Pa, Vietnam 

Lao Cai, Vietnam 

Maolan National Nature Reserve, Libo County, Guizhou, China 
Maolan National Nature Reserve, Libo County, Guizhou, China 
Wuchuan County, Guizhou, China 

Cao Bang , Vietnam 


HNNUDODSIII9 
HNNUO003 
HNNU0607032 
AMNH A161169 
NHMG1401035 
No. AM04005 
KUHE:24635 


KUHE:12898 
HNNUTW9 
KUHE:12896 
FMNH 265919 
FMNH 263415 
HNNUIOOSIO99 
AMNH A163935 
HNNU0606123 
IABHU 21004 
ROM 7472 
ROM 7446 
GZNU20170725018 
CIBjs20171014001 
HNNUOGOSII49 
CIB20130531 
ROM 13044 
HNNUHKOOI 
SYS a002262 
HNNU0606105 
HNNUIO071061 
HNNUIO07D02 
GZNU20170822001 
CIBFJS20150502002 
HNNUO607001 
HNNU707071 
HNNUIOOSIIOI6 
HNNUOO2JL 
HNNUOOIYN 
HNNUOOSIII8S 
HNNU20050032 
20070711017 
ROM 19053 


ROM 38605 paratype 


GZNUJ20180608018 
GZNU20180608019 
HNNUOI9L 
IEBR 3725 


Kuankuoshui National Nature Reserve, Suiyang County, Guizhou, China GZNU2018060802 
Kuankuoshui National Nature Reserve, Suiyang County, Guizhou, China GZNU20180608011 
Kuankuoshui National Nature Reserve, Suiyang County, Guizhou, China GZNU20180608012 


Yuntai Mountain, Shibing County, Guizhou, China 

Nanling Nature Reserve, Ruyuan County, Guangdong, China 
Nanling Nature Reserve, Ruyuan County, Guangdong, China 
Leigongshan Nature Reserve, Leishan County, Guizhou, China 
Longsheng County, Guangxi, China 

Anlong County, Guizhou, China 

Lai Chau, Vietnam 

Ha Giang, Vietnam 

Amami Island, Japan 

Xe Kong, Laos 

Lipu County, Guangxi, China 

Lipu County, Guangxi, China 

Lung Tung Village, Ha Lang , Cao Bang , Vietnam 

Coong Village, Ha Lang , Cao Bang, Vietnam 

Maolan National Nature Reserve, Libo County, Guizhou, China 
Maolan National Nature Reserve, Libo County, Guizhou, China 
Maolan National Nature Reserve, Libo County, Guizhou, China 
Shimian County, Sichuan, China 

Xiling Snow Mountain, Dayi County, Sichuan, China 
Wawushan Mountain, Sichuan, China 

Gai Lai, Vietnam 

Wugongshan Mountain, Jiangxi, China 


Wuyishan, Fujian , China 
Fuzhou City, Fujian, China 
Wuzhishan City, Hainan, China 
Huangshan, Anhui, China 
Hongya, Sichuan, China 


Emeishan City, Sichuan, China 


GZNU20170718002 
CIBHN201108149 
HNNUIOOSIO75 
CIBLS20140616004 
HNNU70028 
HNNUIOOSII09 
AMNH A161439 
AMNH 163782 
IABHU 5275 
FMNH 258107 
NHMGI303018 
NHMG1303019 
IEBR:A2015. 63 
IEBR:A2015. 65 
GZNU20180608007 
GZNU20180608009 
GZNU20160802003 
SM-ZDTW-01 


20130258 
AMI3988 


DYTW-WYS-001 
SCUM-H002CJ 
HNNUOSIII7 
SCUMQ405192C] 
SCUM045199C] 
A-A-WZO001 
HNNU20060103 


KFI85017 
KFI85019 
KF185020 
DQ283345 
MH665669 
DQ835616 
AB200923 
AB511287 
AB200926 
KF185010 
AB200928 


KF185053 
KF185055 
KF185056 
DQ283345 
MH665675 
DQ835616 
AB200947 
AB511287 
AB200950 
KF185046 
AB200952 
DQ650564 
DQ650613 
KF185036 
DQ283394 
KF185038 
AB511284 
AF206487 
AF206484 
MW481360 
MHI93551 
KF185047 
MHI93555 
AF206480 
KF185044 
KT315375 
KF185032 
KFI85041 
KF185052 
MWASI359 
MHI93565 
KF185059 
KF185040 
KF185051 
KF185058 
KF185057 
KF185034 
KF185035 
KF185050 
AF206482 
EUS61556 
MW481353 
MWA81354 
KF185043 
KU356766 
MW481355 
MW481356 
MW481357 
MW481358 
MH193560 
KF185048 
MH193553 
KF185054 
KF185049 
DQ283372 
EU861546 
AB511282 
EUS61542 
MH665676 
KM388701 
5 LCI55910 
LCI55911 
MWA8I350 


KF771294 
KF185000 
DQ283394 
KF185002 
AB511284 
AF206106 
AF206103 
MW481349 
MHI93539 
KFI185011 
MHI193543 
AF206099 
KF185008 
KT315354 
KF184996 
KF185005 
KF185016 
MW481348 
MHI93532 
KF185023 
KF185004 
KFI85015 
KF185022 
KF185021 
KF184998 
KF184999 
KFI185014 
AF206101 


MwW481342 
MW481343 
KF185007 
KU356762 
MWA81344 
MWA81345 
MwW481346 
MW481347 
MHI193540 
KF185012 
MH193533 
KF185018 
KF185013 
DQ283372 


AB511282 


MH665670 


MW481339 
MWwW481340 MWA81351 
MW481341 MW481352 
NC. 029250 NC. 029250 
NC 024180 NC. 024180 
NC 044901 NC. 044901 
KU840519  KU840606 
NC 023949 NC. 023949 
NC 025591 NC 025591 
KX233864 KX233864 
KX269219 KX269219 
KF185031 — KF185067 
KX269222 KX269222 
KX269216 KX269216 
NC 018771 NC. 018771 
KF185029  KF185065 


DQ835616 
AB600991 
AB600990 


DQ650500 
DQ650546 


MWA81371 
MH193605 


MHI93609 
DQ650505 


MWA481370 
MHI93614 


EU861585 
MWA481364 
MW 481365 


MW481366 
MW481367 
MW481368 
MW481369 
MH193615 


MH193607 


EU861572 
AB511282 
EU861568 


MWA81361 
MW481362 
MW481363 
NC. 029250 
NC 024180 
NC 044901 


NC 023949 
NC 025591 
KX233864 


KX269435 
KX269431 
NC_018771 


79 Rana weiningensis Weining County, Guizhou, China SCUMO0405171 KX269217 _KX269217 _KX269432 
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Luo following Fei et al. (2009) and Li et al. (2018a). A total of 
24 morphological characteristics were measured. These 
measurements were as follows: SVL=snout-vent length 
(distance from the tip of the snout to the posterior edge of the 
vent), HDL=head length (distance from the tip of the snout to the 
articulation of jaw), HDW=maximum head width (greatest width 
between the left and right articulations of jaw), SL=snout length 
(distance from the tip of the snout to the anterior corner of the eye), 
IND-internasal distance (minimum distance between the inner 
margins of the external nares), NED=nasal to eye distance 
(distance between the nasal and the anterior corner of 
the eye), NSD=nasal to snout distance (distance between 
the nasal the posterior edge of the vent), IOD-interorbital 
distance (minimum distance between the inner edges of the 
upper eyelids), IFE-distance between anterior corner of 
eye, IAE-distance between posterior corner of eyes, 
UEW-upper eyelid width (greatest width of the upper eyelid 
margins measured perpendicular to the anterior-posterior axis), 
ED-eye diameter (distance from the anterior corner to the posterior 
corner of the eye), TD=tympanum diameter, LAHL=length of 
lower arm and hand (distance from the elbow to the distal end 
of the finger IV), LAN=width of lower arm, HAL-hand length 
(distance from the posterior end of the inner metacarpal tubercle 
to the distal tip of Finger IV), HLL-hindlimb length (maximum 
length from the vent to the distal tip of the toe IV), TL-tibia 
length (distance from knee to tarsus), TW maximal tibia width, 
TFL-length of foot and tarsus (distance from the tibiotarsal 
articulation to the distal end of the toe IV), FOL-foot length (from 
the base of inner metatarsal tubercle to the tip of fourth toe), 
TED=tympanum-eye distance (from anterior edge of 
tympanum to posterior corner of the eye), FDW-finger 
disk width (width at the widest part of the disk of finger 
III, and DPW distal phalanx width (maximal width of 
the distal phalanx of finger III). 

To reduce the impact of allometry, a size-corrected value 
from the ratio of each character to SVL was calculated for 
the following morphometric analyses. Principal component 
analysis (PCAs) of size-corrected variables and simple bivariate 
scatterplots was used to explore and reflect the morphometric 
differences between the new species and O. lipuensis. Mann- 
Whitney U tests were used to test the significance of differences 
on morphometric characters between the new species and O. 
lipuensis, O. kweichowensis and O. wuchuanensis. Mann-Whitney U 
tests were also conducted to test the morphometric differences 
between the males and the females of the new species. The 
statistical analyses were performed using SPSS 21.0 (SPSS, 
Inc, Chicago. IL, USA), and differences were considered to be 
significant at P « 0.05. 

Sex was determined by direct observation of calling behavior 
and the presence of internal vocal sac openings for males, as 
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well as the presence of eggs on the abdomen for females. The 
presence or absence of nuptial pads/spines was examined by 
optical microscopy. 

We compared the morphological characters of the new 
taxon with other species of Odorrana. Comparative data were 
obtained from the literature for 61 species of Odorrana (Table 
3). For comparison, we examined the type and/or topotype 
materials for O. lipuensis, O. kweichowensis, and O. wuchuanensis 
(Appendix 1). 


3. Results 


3.1. Phylogenetic analyses and genetic divergence The 
maximum likelihood (ML) and Bayesian inference (BI) 
phylogenetic trees were constructed based on concatenated 
DNA sequences of the mitochondrial 128 rRNA (752 bp), 16S 
rRNA (1006 bp), and ND2 (1030 bp) genes with a total length 
of 2788 bp. ML and BI analyses resulted in a largely identical 
topology (Figure 2). The new taxon is a sister taxon to O. 
lipuensis with high node support values (100 in BI and 10096 in 
ML; Figure 2). 

The smallest p-distance divergences between the new lineage 
and other species of Odorrana were 6.06% in 12S rRNA (between 
new taxon and O. lipuensis) and 3.90% in 16S rRNA (between 
new taxon and O. geminata), which were at the same divergence 
level as those among recognized congeners (1.94% and 2.79% in 
12S rRNA between O. wuchuanensis vs. O. mutschmanni and O. 
kweichowensis vs. O. schmackeri, respectively; and 130% and 195% 
in 16S rRNA between O. lungshengensis vs. O. yizhanggensis and 
O. nanjiangensis vs. O. hejiangensis, respectively), indicating that 
the new taxon represents an independent evolutionary lineage 
(Tables S2-S3). 


3.2. Morphological analyses The Mann-Whitney U tests 
indicated that males of the new species were significantly 
different from O. lipuensis and O. kweichowensis in many 
morphometric characters (all P values « 0.05; Table 4). PCA 
extracted three and four principal component factors with 
Eigenvalues greater than two in males and females, respectively 
(Table S4). The first two principal components explained 
613496 and 632896, of the total variation in males and females, 
respectively. These differences were mainly influenced by 
limb and head characteristics (Table S4). The new taxon was 
distinctly separated from O. lipuensis and O. kweichowensis on the 
two-dimensional plots of PCI vs. PC2, in both males and females 
(Figure 3). 


3.3. Taxonomic account 

Odorrana liboensis sp. nov. (Tables 2 and 4, Figures 4-5) 
Holotype. GZNU20180608004 (Figure 4), adult male, collected 
by Tao LUO on June 8, 2018 from Maolan National Nature 
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Figure 2 Phylogenetic tree based on mitochondrial 12S rRNA+16S rRNA+ND2 genes. In the phylogenetic tree, ultrafast bootstrap supports (UFB) from 
ML analyses and Bayesian posterior probabilities (BPP) from BI analyses were noted beside nodes. The scale bar represents 0.05 nucleotide substitutions 


per site. Numbers at the tips of branches correspond to the ID numbers in Table 1. 
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Figure 3 Plots of the first principal component (PC 1) versus the second (PC 2) for Odorrana liboensis sp. nov. and O. lipuensis from a principal component 


analysis. A, male. B, female. 


Reserve (25.481711° N, 108.078003" E, ca. 715 m a.s.l.), Libo County, 
Guizhou Province, China. 

Paratypes. Thirteen individuals collected at the same locality 
as the holotype. GZNU20180608002, GZNU20180608003, 
GZNU20180608006, GZNU20180608007, GZNU20180608009, 
GZNU20180608010, GZNU20180608011 and GZNU20180608012, 
adult males collected by Tao Luo on June 8; GZNU20180815001 
(Figure 5), adult females collected by Xiang ZENG on 
August 15, 2018; GZNU20160802001, GZNU20160802002, 
GZNU20160802003, and GZNU20160729009, adult females 
collected by Tao Luo on August 2, 2016. 

Etymology. The specific epithet “liboensis” is in reference to the 
type locality of the new species: Libo County, Guizhou Province, 
China. We propose the common English name “Libo Odorous 
Frog" and Chinese name “Li Bo Chou Wa ( 25i 5E )". 


Differential diagnosis. Odorrana liboensis sp. nov. is assigned to 


genus Odorrana based upon molecular phylogenetic analyses 
and the following morphological characters: (1) dorsum is 
mostly green; (2) tips of digits dilated, tapering, disks with 
circum-marginal grooves or lateroventral grooves, and vertical 
diameter longer than horizontal diameter in the disks; (3) 
supernumerary tubercle below the base of fingers III and IV; 
(4) feet fully webbed to disks, without tarsal fold; (5) the first 
finger thick and nuptial pad distinct and (6) dorsal skin smooth, 
dorsolateral folds absent or fine (Fei et al., 2009; Li et al, 20182). 
Odorrana liboensis sp. nov. is distinguished from congeners by 
a combination of the following characters: (1) having medium 
body size, with the snout-vent length (SVL) of adult females 
approximately 1.2 times as long as that of males at 56.9 + 1.0 


(558—582 mm, n - 9) in females and 487 £ 12 (471-499 mm, 
n = 5) in males; (2) head length greater than width in males 
and females; (3) tympanum distinctly visible, greater than 
one-half the diameter of the eye; (4) eyes big and prominent, 
with the width of the upper eyelid (UEW) approximately 3/4 
of the eye diameter (ED); (5) dorsolateral folds absent; (6) two 
metacarpal tubercles; (7) relative finger lengths: II « I « IV « 
III; (8) subarticular tubercles on fingers prominent: 1, 1, 2, 2; (9) 
one metatarsal tubercle; (10) tibiotarsal articulation reaching 
to between the eye and the nostril when the leg is stretched 
forward; (11) toes with entire webbing to disks; (12) subarticular 
tubercles on toes prominent: 1, 1, 2, 3, 2; (13) dorsal surfaces of 
limbs with distinct brownish-black bands; (14) smooth, grass- 
green dorsum with irregular brown mottling; (15) venter 
smooth, lacking black spots; and (16) lacking pectoral spinules, 
lacking vocal sacs, and a light white nuptial pad present on 
finger I in males. 

Description of holotype. GZNU 20180608004 (Figure 4), 
adult male. Medium body size, SVL 49.8 mm; habitus slender; 
head length slightly larger than head width (HDL=17.6 mm; 
HDW-15.6 mm; HDL/HDW 1.12); snout short, rounded, and 
projecting beyond the lower jaw in dorsal view, longer than 
the diameter of the eye (SL/ED 1.30); nostril rounded, distinct, 
closer to the snout tip than the eye (NEL/NSL 0.52); internasal 
distance less than interorbital distance (IND/IOD 0.91) and 
greater than upper eyelid width (IND/UEW 1.06); pineal body 
invisible; tympanum distinct, rounded, 7596 of eye diameter, 
depressed relative to the skin of the temporal region, tympanic 
rim slightly elevated relative to tympanum; vomerine teeth 
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Table 3 References for morphological characters for congeners of the genus Odorrana. 


Ka 


Species 


Literature obtained 


Odorrana absita (Stuart & Chan-ard, 2005) 

Odorrana amamiensis (Matsui, 1994) 

Odorrana andersonii (Boulenger, 1882) 

Odorrana anlungensis (Liu & Hu, 1973) 

Odorrana aureola Stuart, Chuaynkern, Chan-ard & Inger, 2006 


Odorrana bolavensis (Stuart & Bain, 2005) 
Odorrana cangyuanensis (Yang, 2008) 

10 Odorrana chapaensis (Bourret, 1937) 

11 Odorrana chloronota (Günther, 1876) 

12 Odorrana dulongensis Liu, Che & Yuan, 2021 
13 Odorrana exiliversabilis Li, Ye & Fei, 2001 


NO 0o M Ov Un 4 Q9 t2 eC 


14 Odorrana fengkaiensis Wang, Lau, Yang, Chen, Liu, Pang & Liu, 2015 


15 Odorrana geminata Bain, Stuart, Nguyen, Che & Rao, 2009 
16 Odorrana gigatympana (Orlov, Ananjeva & Ho, 2006) 

17 Odorrana grahami (Boulenger, 1917) 

18 Odorrana graminea (Boulenger, 1900) 

19 Odorrana hainanensis Fei, Ye & Li, 2001 

20 Odorrana hejiangensis (Deng & Yu, 1992) 

21 Odorrana hosii (Boulenger, 1891) 

22 Odorrana huanggangensis Chen, Zhou & Zheng, 2010 

23 Odorrana ichangensis Chen, 2020 

24 Odorrana indeprensa (Bain & Stuart, 2006) 

25 Odorrana ishikawae (Stejneger, 1901) 

26 Odorrana jingdongensis Fei, Ye & Li, 2001 

27 Odorrana junlianensis Huang, Fei & Ye, 2001 

28 Odorrana khalam (Stuart, Orlov & Chan-ard, 2005) 

29 Odorrana kuangwuensis (Liu & Hu, 1966) 

30 Odorrana kweichowensis Li, Xu, Lv, Jiang, Wei & Wang, 2018 
31 Odorrana leporipes (Werner, 1930) 

32 Odorrana lipuensis Mo, Chen, Wu, Zhang & Zhou, 2015 

33 Odorrana livida (Blyth, 1856) 

34 Odorrana lungshengensis (Liu & Hu, 1962) 

35 Odorrana macrotympana (Yang, 2008) 

36 Odorrana margaretae (Liu, 1950) 

37 Odorrana mawphlangensis (Pillai & Chanda, 1977) 

38 Odorrana monjerai (Matsui & Jaafar, 2006) 

39 Odorrana morafkai (Bain, Lathrop, Murphy, Orlov & Ho, 2003) 


40 Odorrana mutschmanni Pham, Nguyen, Le, Bonkowski & Ziegler, 2016 


41 Odorrana nanjiangensis Fei, Ye, Xie & Jiang, 2007 
42 Odorrana narina (Stejneger, 1901) 

43 Odorrana nasica (Boulenger, 1903) 

44 Odorrana nasuta Li, Ye & Fei, 2001 

45 Odorrana orba (Stuart & Bain, 2005) 

46 Odorrana rotodora (Yang & Rao, 2008) 

47 Odorrana schmackeri (Boettger, 1892) 

48 Odorrana sinica (Ahl, 1927) 


49 Odorrana splendida Kuramoto, Satou, Oumi, Kurabayashi & Sumida, 2011 


50 Odorrana supranarina (Matsui, 1994) 

51 Odorrana swinhoana (Boulenger, 1903) 

52 Odorrana tianmuii Chen, Zhou & Zheng, 2010 
53 Odorrana tiannanensis (Yang & Li, 1980) 

54 Odorrana tormota (Wu, 1977) 

55 Odorrana trankieni (Orlov, Le & Ho, 2003) 

56 Odorrana utsunomiyaorum (Matsui, 1994) 

57 Odorrana versabilis (Liu & Hu, 1962) 

58 Odorrana wuchuanensis (Xu, 1983) 

59 Odorrana yentuensis Tran, Orlov & Nguyen, 2008 
60 Odorrana yizhangensis Fei, Ye & Jiang, 2007 
61 Odorrana zhaoi Li, Lu & Rao, 2008 


moderately developed, on two oblique ridges; tongue cordiform, 
deeply notched posteriorly; and eyes large, slightly protuberant 
in dorsal view, eye diameter 26.70% of head length, pupils 
transverse; supratympanic fold absent. 

Forelimbs are slender and comparatively short, the length 
of the lower arm and hand is 419796 of SVL; fingers slender, 
relative finger lengths: I«IIXIV«IIT tips of all except first finger 


Odorrana bacboensis (Bain, Lathrop, Murphy, Orlov & Ho, 2003) 
Odorrana banaorum (Bain, Lathrop, Murphy, Orlov & Ho, 2003) 


Stuart and Chan-ard (2005) 
Matsui (1994) 
Boulenger (1882) 

Hu et al. (1973) 

Stuart et al. (2006) 
Bain et al. (2003) 

Bain et al. (2003) 
Stuart and Bain (2005) 
Yang (2008) 

Bourret (1937) 
Günther (1876) 

Liu et al. (2021) 

Fei et al. (2001b) 
Wang et al. (2015) 
Bain et al. (2009) 
Orlov et al. (2006) 
Boulenger (1917) 
Boulenger (1900) 

Fei et al. (20012) 

Deng et al. (1992) 
Boulenger (1891) 
Chen et al. (20102) 
Shen et al. (2020) 
Bain and Stuart (2005) 
Stejneger (1901) 

Fei et al. (20012) 

Ye and Fei (2001) 
Stuart and Chan-ard (2005) 
Hu et al. (1966) 

Li et al. (20182) 
Werner (1930) 

Mo et al. (2015) 

Blyth (1856) 

Liu and Hu (1962) 
Yang (2008) 

Liu (1950) 

Pillai Chanda (1977) 
Matsui and Jaafar (2006) 
Bain et al. (2003) 
Pham et al. (2016) 

Fei et al. (20072) 
Stejneger (1901) 
Boulenger (1903) 

Fei et al. (2001b) 
Stuart and Bain (2005) 
Yang (2008) 

Boettger (1892) 

Ahl (1927) 

Kuramoto et al. (2011) 
Matsui (1994) 
Boulenger (1903) 
Chen et al. (2010b) 
Yang and Li (1980) 
Wu (1977) 

Orlov et al. (2003) 
Matsui (1994) 

Li and Hu (1962) 

Wu et al. (1983) 

Tran et al. (2008) 

Fei et al. (2007b) 

Li et al. (2008) 


expanded with circummarginal grooves, horizontal grooves 
present, without lateral fringes; width of finger III disc about 
68.57% the diameter of tympanum; webbing absent; subarticular 
tubercles prominent: 1, 1, 2, 2; inner metacarpal tubercle oval, 
elongate; outer metacarpal tubercle oval; and light white nuptial 
pads present. 

Hindlimbs are slender (HLL/SVL 1.86); heels overlapping 
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when thighs are positioned at right angles to the body; 
tibiotarsal articulation reaching the between eye to nostril 
when leg stretched forward; foot length less than tibia length 
(FL/TL 0.87); relative toe lengths I<II<III<V<IV; tips of toes 


expanded into disc with circummarginal grooves; subarticular 
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tubercles prominent: 1, 1, 2, 3, 2; toes with entire webbing to 
disks; elongate, oval inner metatarsal tubercle, larger than toe I; 
and outer metatarsal tubercle absent. 

Skin is smooth on upper surfaces; venter smooth; tiny 
spinules on lateral body, upper edge of lid, temporal region, and 


Figure 4 Morphological features of the live adult male holotype GZNU20180608004 of Odorrana liboensis sp. nov. and O. wuchuanensis. (A) Dorsal view; 
(B) Ventral view; (C) Dorsolateral view; (D) Tadpole; (E) Type locality of Odorrana liboensis sp. nov. in June, 2018 (dry season); (F) Dorsolateral view of O. 


wuchuanensis in life. 


INo.4 


anterior and posterior edge of tympanum; weak supratympanic 
fold from the posterior edge of the eye to the posterior edge of 
the tympanum; and dorsolateral fold absent. 

Coloration of holotype in life (Figure 4). Dorsal parts of 
head and dorsum, flank, forelimb, thigh, tibia, and foot are 
grass-green with irregular brown mottling; the throat, chest, 
and abdomen lack black spots; dorsal surfaces of limbs have 
indistinct brownish-black bands, the thighs with four brownish- 
black bands and tibias with three; ventral surface of the limbs is 
pink; the iris is black, surrounded by a gold-green network; and 
the tympanum is dark brown. 
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Preserved holotype coloration. After preservation in 7596 
ethanol, the dorsal surface of the body coloration changed to 
dark brown grey; the dorsal surface of the head coloration 
changed to dark grey; the transverse bands on limbs and digits 
were not distinct and the coloration changed to lighter colors; 
the throat was light grayish yellow; the chest was light ash 
black; the belly was light gray-white; the posterior of ventral 
surface of body, inner thigh, and upper part of the tibia were 
creamy yellow; the palms and metatarsal tubercles were ash 
black; and the ventral surface of the forelimbs and hindlimbs 


were creamy yellow with brown mottling. 


Table 4 Morphological comparison of Odorrana liboensis sp. nov. (OB), O. lipuensis (OP), O. kweichowensis (OK), and O. wuchuanensis (OW). All units are in 
mm. P-values are at 95% significance. Morphometric characters are explained in the methods section. BM and BF are the abbreviations for male and female 


from Odorrana liboensis sp. nov. 


Odorrana liboensis sp. nov. 


Odorrana lipuensis 


Odorrana wuchuanensis 


Measure- 


m Male (n = 5) Female (n =9) Male (n =4) Female (n =4) Male (n =12) Female (n =5) 
Range Mean+SD Range Mean+SD Range MeantSD Range Mean+SD Range Mean+SD Range Mean+SD 
SVL 471-499 487412 558-582 569410 407-477 431431 499-541 522418 700-805 1019+25 998-1057 768+32 
HDL 167-184 175406 193-215 207408 146-156 150404 172-197 186412 226-289 324416 299-341 262417 
HDW 138-156 149406 182-195 189405 132-146 137406 173-205 183415 217-282 338438 309-400 244116 
SL 65-82 73405 85-95 90x04 566-66 61405 79-82 81x02 94-178 132419 103-155 126422 
IND 48-56 52:03 53-68 61407 36-45 42x04 45-55 50:05 74-100 105405 98-111 87 £07 
NED 31-45 37 +05 52-66 5905 35-43 3.7 +04 46-52 48 +03 76-85 91x05 83-96 80x03 
NSD 28-4 34404 27-33 31403 22-32 28 +05 29-37 33 +04 39-73 5.3 +06 43-59 47 +10 
IOD 46-68 60+08 72-79 75 +03 32-41 3.7 +04 47-56 50+04 58-135 6.6412 51-83 802.1 
IFE 83-9.6 93+04 102-122 112 +08 68-83 Z507 85-95 90:04 124-159 1146+11 105-133 136411 
IAE 143-153 147403 152-173 160409 99-114 105406 129-149 140+09 225-299 294+08 285-305 269427 
UEW 35-48 41:04 44-51 48 +03 32-40 3.6 +03 47-52 50+02 45-105 73E l8 43-90 64415 
ED 47-65 5.6 + 0.6 55-69 64405 43-55 47 +06 51-65 59106 67-133 9615 69-107 86417 
TD 31-38 35x03 48-59 51x04 32-41 3.6 +04 38-51 44x06 46109 65+19 35-84 64416 
LAHL 205-226 212406 282-297 290407 178-224 197421 256-278 265£09 323-417 491424 464-518 386425 
LAN 3-47 37 £06 35-5] 42x06 37-50 43:06 32-48 37:08 46-109 91405 84-97 83+16 
HAL 12-128 124405 121-144 132+10 103-128  116X10 141-158 147408 248-295 294405 289-299 279+15 
HLL 868-1027 927459 1027-1115 106840 682-769 731437 262-279 27A+08 1035-1402 1647+75 1589-1777 1290+388 
TL 264-358 307432 311-415 363445 216-245 233413 265-289 277410 335-482 517407 507-525 446437 
TW 63-69 66x02 63-76 70x05 39-56 48 +07 46-66 53+09 78-424 93403 89-98  22.0+168 
TFL 354-387 373412 422-485 442426 309-338 318413 261-386 350460 403-487 444420 409-459 447425 
FOL 245-269 262407 272-313 289417 195-230 213415 254-268 263406 324-463 SI5+20 485-534 421439 
TED 19-22 21401 21-34 2705 13-15 14+01 18-23 21x02 29-46 45 +08 34-51 37 +05 
FDW 21-25 23202 25-31 28:03 12-20 16303 17-26 21404 31-67 49x05 42-55 44211 
DPW 25-28 2601 22-28 24t03 11-16 14402 17-23 2.0 +03 14-77 59+05 52-65 5116 
HDL/SVL 0.34-0.39 036+0.02 0.33-0.37 036+0.02 033-036 035+002 033-037 036+0.02 032-038 032+001 030-033 0.34 0.02 
HDW/SVL 028-033 0314002 032-034 0334001 031-033 0324001 033-038 0.35+0.02 030-037 033+0.03 031-039 0324002 
HDL/HDW 113-123 117x004 101-115 110+0.05 107-114 110+0.03 096-107 102+005 099-115 0.97+008 085-106 108+0.05 
SL/SVL 013-017 015x001 015-016 016+000 013-016 0.14+001 015-016 016+0.00 013-024 0134002 010-015 0.16 +0.03 
SL/HDL 038-045 042003 041-049 04345003 038-045 0414003 042-048 0443003 042-061 040-004 034-045 0.48 +0.06 
IOD/HDW 031-047 040-006 039-041 0404001 024-031 027:003 023-031 0274003 025-048 019x003 016-022 032:t007 
IND/IOD 073-121 089x015 074-092 0814007 088-13] 113x018 084-115 100x015 057-142 1634029 133-210 114x023 
ED/TD 134-202 1634021 110-135 126x013 123-150 1334012 116-144 135x013 121-152 1534029 122-198 135011 
ED/HDL 027-036 0323003 029-032 0314001 028-035 03145003 030-034 03242002 026-046 0293x004 023-033 0332005 
TD/HDL 017-022 020x001 023-028 025x002 022-026 0244002 022-026 024-002 017-038 0204005 012-025 0244x005 
HAL/SVL 0.22-0.27 025x001 022-026 0234002 025-029 0274002 027-029 0284001 033-041 029x000 028-029 036:003 
HLL/SVL 175-217 1914014 183-192 1884005 156-185 170x012 051-053 052x001 128-185 1624009 154-178 1684014 
TED/HDL 011-013 012x001 010-016 0134002 008-010 009x001 010-012 011x001 011-017 0144003 010-016 0144002 
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(Continued Table 4) 
Odorrana kweichowensis P-value from Mann-Whitney U test 
id Male (n = 18) Female (n = 17) Male Female 
Range  MeantSD Range  MeantSD OBvs.OP OBvs.OW OBvs.OK  OBwOP OBvs.OW OBvs.OK BM 1s. BF 
SVL 424-488  448+18 794-940 855447 0.007 0,000 0.000 0016 0.008 0.000 0.001 
HDL 151-179 16.3 +08 257-314 280 + 16 0.001 0.034 0.463 0.111 0.008 0.002 1.000 
HDW 113-134 120205 207-251 DET 0.007 0382 0.000 0.063 0421 0,000 0004 
SL 6.0-7.7 70+05 117-160 128 +11 0.004 0.193 0.131 0.730 0.008 0.101 0112 
IND 42-57 50x05 81-118 SE 0.001 0.129 0.322 0.190 0.548 0.820 1000 
NED 32-43 38 +03 61-93 72+09 1000 0.000 0.053 0.111 0.032 0.001 0.001 
NSD 31-39 33402 53-78 59 +06 0.019 0.041 0,059 0.016 0.690 0.000 0.001 
IOD 39-54 44t04 39-107 813:15 0.001 0.015 0.000 0.016 0.008 0,000 0364 
IFE 40-77 59 t10 99-153 POTIS 0.001 0.082 0.000 0.016 0.008 0.000 0.052 
IAE 5.6-10.0 79216 165-208 18.2 +13 0.001 0.001 0.000 0.413 0.421 0,000 0.042 
UEW 35-46 39403 45-59 53404 0.190 0.422 0527 0016 0.056 0.000 0.797 
ED 49-65 5.5 +05 81-102 90x05 0.042 0310 0145 1000 0.095 0.048 0797 
TD 31-44 36403 38-54 46x04 0.699 0.018 0.020 0.286 0.008 0.000 0.001 
LAHL 204-250 22.4 +12 386-450 415 +18 0.298 0.001 0,000 0.905 0.095 0.120 0.001 
LAN 33-56 45+06 62-96 72x09 0.019 0.001 0001 0.556 0,032 0.048 0.898 
HAL 112-145 131x 09 211-260 238 £13 0.240 0.000 0,000 0.016 0,008 0,000 0.029 
HLL 752-875 816436 1393-1611 1523466 0.001 0.000 0.194 0.016 0.008 0.015 0.797 
TL 217-291 247 £18 472-540 50.0 +19 0.001 0.129 0.002 0.032 0,008 0.319 1.000 
TW 44-65 51x05 93-132 10512 0.001 0554 0.000 0.111 0.008 1.000 0.019 
TFL 32.1-38.9 357 £19 635-730 683 + 28 0.001 0,000 0.085 0.016 0,008 0.085 1000 
FOL 216-263 238+14 . 380-477 432426 0.001 0,082 0.232 0.730 0841 0.880 0.083 
TED 17-25 20+ 0.2 28-39 34403 0.001 0.095 0.212 0.286 0.690 0.085 0.797 
FDW 13-21 17x02 24-36 29 +03 0.001 0.111 0.000 0.190 0.841 0.000 0.699 
DPW 13-19 16+02 22-33 2, X05 0,001 0.001 0,000 0.190 0,008 0.000 0.001 
HDL/SVL 033-038 0364001 031-035 033+ 0.01 0.898 0.034 0.463 0.111 0.008 0.002 1.000 
HDW/SVL 0.25-0.29 0.27+001 024-027 0.26 +0.01 0.083 0.382 0.000 0.063 0.421 0.000 0.004 
HDL/HDW 122-147 136+006 117-130 125x003 0.007 0.000 0.000 0.063 0.032 0.000 0.007 
SL/SVL 0.12-0.18 0.16 + 0.01 014-017 015 € 0.01 0.699 0.193 0.131 0.730 0,008 0.101 0.112 
SL/HDL 034-049 0434004 044-051 0.46 +0.02 0.898 0.002 0253 1000 0421 0.039 0518 
IOD/HDW . 033-040 037-002 019-047 036-006 0,001 0.015 0.085 0.016 0.008 0.410 0.364 
IND/IOD 091-143  l1l4t015 091-217 119+028 0.029 0.006 0.000 0.063 0.008 0.000 0.438 
ED/TD 131-182 155 + 0.18 158-254 196 + 020 00D 0.001 0322 0.413 0.056 0.000 0.007 
ED/HDL 029-039 034x003 028-035 0.32 + 0.02 1000 0972 0.118 0.730 0548 0.140 0.438 
TD/HDL 018-025 0224002 013-019 017 + 001 0.002 0.001 0.005 0283 0.056 0.000 0.001 
HAL/SVL 026-034 0294002 026-033 0.28 + 0.02 0.060 0.000 0,000 0.016 0.008 0.000 0.029 
HLL/SVL 163-202 182 + 010 167-197 178 + 008 0.001 0.000 0194 0.016 0.008 0.015 0797 
TED/HDL 010-014 013x001 010-014 0.12 +0.01 0.001 0.003 0176 0286 0.690 0.543 0.438 


Variations. The basic statistics for measurements are presented 
in Table 2. In life, all paratypes matched the overall basic 
morphological characters of the holotype. Females (SVL 569 + 
10 mm, n = 5) have larger body size than males (SVL 487 + 12 
mm, n = 9), with the SVL in females approximately 1.2 times 
that in males, but the relative sizes of ED and TD are obviously 
larger in males than in females (Table 4); the dorsal surfaces of 
limbs with indistinct brownish-black bands (Figure 5). 

Sexual dimorphism. Adult females have an SVL of 558-582 
mm, larger than adult males, with an SVL of 47.1-499 mm. 
Adult males lack vocal sacs. In breeding, light white nuptial pads 
are present on finger I in males. 


Comparisons. Comparative data of Odorrana liboensis sp. nov. 
with 61 recognized congeners of Odorrana are given in Table S5. 
By having medium body size (minimum SVL » 47.0 mm 
in males), Odorrana liboensis sp. nov. differs from O. absita, O. 
anlungensis, O. gigatympana, O. huanggangensis, O. khalam, 
O. kweichowensis, O. monjerai, O. morafkai, O. nasica, O. orba, 
and O. tormota (vs. maximum SVL « 47.0 mm in males). By 
having medium body size (maximum SVL « 50.0 mm in 
males), Odorrana liboensis sp. nov. differs from O. amamiensis, O. 
andersonii, O. aureola, O. bacboensis, O. cang yuanensis, O. chapaensis, 
O. geminata, O. grahami, O. hainanensis, O. hosii, O. ishikawae, 
O. jingdongensis, O. junlianensis, O. kuangwuensis, O. leporipes, 
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O. livida, O. lungshengensis, O. macrotympana, O. margaretae, O. 
mutschmanni, O. nanjiangensis, O. nasuta, O. sinica, O. splendida, O. 
supranarina, O. swinhoana, O. tiannanensis, O. trankieni, O. versabilis, 
O. wuchuanensis, and O. zhaoi (vs. minimum SVL > 50.0 mm). 

By having medium body size (vs. maximum SVL < 59.0 mm 
in females), Odorrana liboensis sp. nov. differs from O. amamiensis, 
O. andersonii, O. anlungensis, O. aureola, O. bacboensis, O. banaorum, 
O. bolavensis, O. chapaensis, O. chloronota, O.dulongensis, O. 
fengkaiensis, O. geminata, O. grahami, O. hainanensis, O. hejiangensis, 
O. hosii, O. huanggangensis, O. ishikawae, O. jingdongensis, O. 
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junlianensis, O. kuangwuensis, O. kweichowensis, O. livida, O. 
lungshengensis, O. macrot ympana, O. margaretae, O. monjerai, O. 
mora f kai, O. mutschmanni, O. nanjiangensis, O. narina, O. nasuta, 
O. orba, O. rotodora, O. schmackeri, O. splendida, O. supranarina, O. 
swinhoana, O. tianmuii, O. tiannanensis, O. tormota, O. versabilis, O. 
wuchuanensis, O. yentuensi, O. yizhangensis, and O. ichangensis (vs. 
maximum SVL > 59.0 mm). 

By lacking dorsolateral folds, Odorrana liboensis sp. nov. 
differs from O. absita, O. amamiensis, O. banaorum, O. bolavensis, O. 


exiliversabilis, O. gigat ympana, O. graminea, O. hosii, O. indeprensa, 


Figure 5 Morphological features of the live adult female GZNU20180815001 of Odorrana liboensis sp. nov. (A) Dorsal view; (B) Dorsolateral view; (C) 
Ventral view; (D) Ventral view of right foot; (E) Ventral view of right hand; (F) The entrance habitat of the karst cave at the type locality of Odorrana 


liboensis sp. nov.; (G) Water outlet in the cave. 
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O. khalam, O. leporipes, O. livida, O. monjerai, O. narina, O. nasica, 
O. nasuta, O. orba, O. supranarina, O. tormota, O. trankieni, O. 
utsunomi yaorum, O. yentuensis, and O. zhaoi (vs. present). 

By the tibiotarsal articulation reaching to between the eye 
and the nostril when the leg is stretched forward, Odorrana 
liboensis sp. nov. differs from O. bacboensis, O. jingdongensis, O. 
lungshengensis, O. margaretae, O. mutschmanni, O. nanjiangensis, 
O. narina, O. orba, O. sinica, O. swinhoana, O. tormota, and O. 
yizhangensis (vs. reaching the tip of the snout), from O. nasica and 
O. nasuta (vs. reaching the tip of the snout or a little beyond), 
from O. hainanensis (vs. reaching the tip of the snout or the 
anterior corner of eye), from O. junlianensis (vs. reaching the 
tip of the snout or between the nostril and the snout), from O. 
cang yuanensis, O. exiliversabilis, O. fengkaiensis, O. gigatympana, 
O. grahami, O. graminea, O. tiannanensis, O. versabilis, O. yentuensis, 
and O. zhaoi (vs. reaching to or beyond the tip of the snout), from 
O. amamiensis (vs. reaching far beyond the tip of the snout), from 
O. amamiensis, O. anlungensis, O. huanggangensis, O. kuangwuensis, 
O. macrot ympana, O. wuchuanensis, and O. ichangensis (vs. reaching 
the nostril or beyond the tip of the snout), from O. lipuensis, O. 
splendida, and O. supranarina (vs. reaching the anterior corner of 
the eye), from O. rotodora and O. trankieni (vs. reaching beyond 
the eye), and from O. utsunomiyaorum (vs. reaching between the 
anterior corner of the eye and the nostril). 

By lacking black bars on the lips, Odorrana liboensis sp. 
nov. differs from O. andersonii, O. anlungensis, O. bacboensis, O. 
bolavensis, O. chapaensis, O. dulongensis, O. fengkaiensis, O. geminata, 
O. grahami, O. hainanensis, O. hejiangensis, O. huanggangensis, O. 
ishikawae, O. junlianensis, O. kuangwuensis, O. kweichowensis, O. 
lungshengensis, O. margaretae, O. mutschmanni, O. nanjiangensis, O. 
schmackeri, O. splendida, O. supranarina, O. tianmuii, O. tiannanensis, 
O. wuchuanensis, O. yizhangensis, and O. ichangensis (vs. the 
presence of black bars). 

By the absence of vocal sacs in males, Odorrana liboensis 
sp. nov. differs from O. absita, O. amamiensis, O. andersonii, O. 
anlungensis, O. aureola, O. bacboensis, O. banaorum, O. bolavensis, 
O. cang yuanensis, O. chapaensis, O. chloronota, O. dulongensis, 
O. exiliversabilis, O. fengkaiensis, O. geminata, O. gigatympana, 
O. grahami, O. graminea, O. hainanensis, O. hejiangensis, O. 
huanggangensis, O. indeprensa, O. ishikawae, O. jingdongensis, O. 
junlianensis, O. khalam, O. kweichowensis, O. lungshengensis, O. 
macrotympana, O. mora f kai, O. nanjiangensis, O. nasica, O. nasuta, 
O. orba, O. sinica, O. swinhoana, O. tianmuii, O. tiannanensis, O. 
tormota, O. trankieni, O. utsunomi yaorum, O. versabilis, O. yentuensis, 
O. yizhangensis, O. ichangensis, and O. zhaoi (vs. the presence of 
vocal sacs). 

By the relative lengths of fingers I<II<IV<III, Odorrana 
liboensis sp. nov. differs from O. lipuensis and O. dulongensis (vs. 
I-II-IV«IIT) from O. absita, O. amamiensis, and O. anlungensis(vs. 
IIcI-IV«IIT) from O. andersonii (vs. I<l-IV<M), from O. aureola, O. 
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bacboensis, O. banaorum, O. bolavensis, O. chloronota, O. f engkaiensis, 
O. gigatympana, O. grahami, O. jingdongensis, O. junlianensis, O. 
kweichowensis, O. margaretae, O. mutschmanni, O. orba, O. schmackeri, 
O. swinhoana, O. tiannanensis, O. wuchuanensis, and O. yentuensis 
(vs. IIXIXIV«IID, from O. cang yuanensis (vs. IXIV«II«IIT), from 
O. exiliversabilis, O. hainanensis, O. hejiangensis, O. kuangwuensis, 
O. nasica, O. nasuta, and O. versabilis (vs. IIXIV«I«IID, from O. 
huanggangensis (vs. Is1I=IV<IIl), from O. splendida, O. supranarina, 
and O. trankieni (vs. IV<II<I<II), from O. zhaoi (vs. I-IV <<). 

By lacking white pectoral spinules in mature males, Odorrana 
liboensis sp. nov. differs from O. andersonii, O. fengkaiensis, O. 
grahami, O. graminea, O. hainanensis, O. huang gangensis, O. 
jingdongensis, O. junlianensis, O. kweichowensis, O. lungshengensis, O. 
margaretae, O. tianmuii, O. yizhangensis, and O. zhaoi (vs. presence 
of white pectoral spinules). 

By having two metacarpal tubercles on the base of the 
hand, Odorrana liboensis sp. nov. differs from O. absita, O. aureola, 
and O. bacboensis (vs. lacking metacarpal tubercles), from O. 
cangyuanensis, O. dulongensis, O. lungshengensis, O. nanjiangensis, 
O. orba, O. schmackeri, O. tiannanensis, and O. yizhangensis (vs. 
one metacarpal tubercle), from O. andersonii, O. exiliversabilis, O. 
fengkaiensis, O. hainanensis, O. hejiangensis, O. huanggangensis, O. 
jingdongensis, O. junlianensis, O. nasica, O. nasuta, O. tianmuii, O. 
tormota, and O. ichangensis (vs. three metacarpal tubercles). 

By having one metacarpal tubercle on the base of the 
metatarsus, Odorrana liboensis sp. nov. differs from O. absita, 
O. amamiensis, O. exiliversabilis, O. nasica, and O. nasuta (vs. two 
metatarsal tubercles). 

The congeners O. kweichowensis and O. wuchuanensis have a 
sympatric distribution with Odorrana liboensis sp. nov. The new 
species can be distinguished from these species by a series of 
morphological characters as follows. The new species differs 
from O. kweichowensis by having a larger body size (adult males 
with a length of 47.1-499 mm vs. adult O. kweichowensis males 
with a length of 362-433 mm) and having a lower ratio of 
IND/IOD at 0.89 in males and 0.81 in females (vs. 114 in males 
and 119 in females for O. kweichowensis), lacks black bars on the 
lips (vs. the presence of black bars), an absence of vocal sacs in 
males (vs. the presence of vocal sacks), lacks large black spots on 
the dorsum (vs. large black spots in the center of the dorsum), 
has thighs with four brown bands and tibias with three (vs. 
thighs with five brown bands and tibias with six), lacks white 
pectoral spinules in mature males (vs. the presence of white 
pectoral spinules), and lacks toes with entire webbing (vs. full 
webbing). 

The new species differs from O. wuchuanensis by having 
a smaller body size (SVL 47.1-499 mm in adult males and 
55.8-582 mm in adult females vs. 71.1-765 mm in adult males 
and 75.8-900 mm in adult females) and a lower ratio of IND/ 
IOD 039 in males and 0.81 in females (vs. 163 in males and 114 
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in females in O. wuchuanensis). The new species has a tibiotarsal 
articulation reaching to between the eye and the nostril when 
the leg is stretched forward (vs. reaching the nostril), lacks black 
bars on the lips (vs. the presence of black bars), lacks large black 
spots on the ventral surface (vs. the presence of large black 
spots), has relative finger lengths of LIN (vs. IT-I-IV«TIT), 
lacks black spots on the dorsum and ventral surfaces (vs. having 
large black spots in the dorsum and ventral surfaces), and 
features males without white spines on the dorsal surface of the 
arm (vs. the presence of large white spines). 

Odorrana liboensis sp. nov. is phylogenetically closest to O. 
lipuensis, and this new species could be distinguished from the 
latter by tibiotarsal articulation reaching to between the eye 
and the nostril when the leg is stretched forward (vs. reaching 
the anterior corner of the eye); relative finger lengths I<II<IV<III 
(vs, I-IIXIV«IID; males and females with a lower ratio of TD/ED 
(mean 0.89 vs. O. lipuensis, mean 1.15); males and females with a 
higher ratio of IOD/HDW (in males 0.31-0.47, mean 0.40, and 
in females 0.39-0.41, mean 0.40 vs. in males 0.24-0.31, mean 027, 
and in females 023-031, mean 027); males with a higher ratio 
of TED/HDL (0.11-0.13, mean 0.12 vs. 0.08-0.10, mean 0.09); the 
absence of small white spiny grains from the anterior corner of 
the eye along the eye and via the tympanic membrane below 
to the anterior part of the cloacal foramen (vs. the presence of 
grains); indistinct brownish-black bands on limbs (vs. distinct); 
the throat, chest, and abdomen lacking black spots (vs. the 
presence of black spots); and the ventral surface of the limbs is 
pink (vs. light purple). 

Distribution and ecology. Odorrana liboensis sp. nov. is known 
only from the type locality, Maolan National Nature Reserve, 
Libo County, Guizhou Province, China, at elevations between 
645 and 728 m. The new species has only been found in one cave 
in the area located far from the village. There is no light in the 
cave, and the new species was found in a small pool approximately 
15 m wide and 79 cm deep, about 145 m from the entrance of 
the cave, where the water temperature is approximately 20 °C all 
year round. The tadpoles were collected on July 23, 2016, but no 
pairs of male and female adults were found to hold them, and 
with four years of survey data, the adults were collected only 
from mid-July to mid-August. Therefore, we speculate that the 
breeding period begins in late June and continues until about mid- 
August. Inside this cave, Odorrana liboensis sp. nov. is sympatric 
with Chinapotamon maolanense, Sinoc yclocheilus longibarbatus, 
Hipposideros armiger, Aselliscus stoliczkanus, Murina liboensis, and 
Leopoldamys edwardsi. Outside the cave, no adults and tadpoles 
of the new species were found during a herpetological survey 
in the vicinity. However, a number of amphibians and reptiles 
can be found outside the cave, including Tylototriton asperrimus, 
Quasipaa boulengeri, Bufo gargarizans, Kurixalus eif fingeri, 
Rhacophorus denn ysi, Goniurosaurus liboensis, Sinomicrurus 
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macclellandi, Lycodon flavozonatum, and Orthriophis moellendor f fi. 
4. Discussion 


Most species of Odorrana live in montane streams. Previously, 
only two species (O. wuchuanensis and O. lipuensis) had been 
reported to live in cave environments (Fei et al., 2012; Mo et a, 
2015). We discovered, and report on, a third species of Odorrana, 
Odorrana liboensis sp. nov., that inhabits cave environments. 
Two of these species (including O. wuchuanensis) are found in 
karst caves in Guizhou. Phylogenetic analyses based on three 
mitochondrial genes suggested that Odorrana liboensis sp. nov. 
belongs to Odorrana but is distinct from its congeners. The 
genetic distances of 128 rRNA and 16S rRNA between the 
Odorrana liboensis sp. nov. and the closely related O. lipuensis were 
6.06% and 5.1996. The genetic distances of 16S rRNA was greater 
than 596, and this is greater than the distance that typically 
represents differentiation, at the species level, in frogs (2396) 
(Vences et al, 2005; Fouquet et al., 2007). The new species was 
different from its congeners on the basis of many morphological 
characters, and this supports its validity. Odorrana liboensis sp. 
nov. described here increases the number of Odorrana species to 
62, with 39 recorded from China (Frost, 2021; AmphibiaChina, 
2021). 

Based on 49 previously named species, one new species, and 
three mitochondrial genes, we conducted a phylogenetic study 
of the genus Odorrana. The number of species covered here and 
the amount of data analyzed exceeds previously reported data 
(Chen et a, 2013; Li et al, 2018a). Eight highly supported major 
clades were identified in Odorrana (Clades A-H). The topology 
of this tree differed significantly from previous studies, focusing 
mainly on the root evolutionary branches of Odorrana (Figure 
2). In the phylogenetic tree, O. lipuensis + Odorrana liboensis sp. 
nov. was the first species to diverge from the genus Odorrana, 
whereas in other studies by Chen et al. (2013) and Li et al. 
(2018a), the first species to diverge was O. chapaensis. In our 
phylogenetic tree, O. chapaensis + O. geminata (Clade D1) serves 
as a sister taxonomic unit to Clade D2 consisting of O. andersonii 
as well as O. wuchuanensis, O. mutschmanni, O. yizhangensis, O. 
lungshengensis, and O. anlungensis. In contrast, in previous studies, 
O. chapaensis appeared as a sister taxon in a clade that included 
all other Odorrana (Chen et al., 2013; Li et al., 2018a), or as a sister 
taxon to Odorrana other than O. lipuensis (He, 2017). Ye and Fei 
(2001) suggested that the primitive taxa of Odorrana may have 
originated from the Hengduan Mountains and the plateau of 
western Yunnan. The Guizhou plateau may have been the center 
of differentiation of Odorrana (Ye and Fei, 2001). O. lipuensis and 
Odorrana liboensis sp. nov. appear to have diverged from Odorrana 
and formed the ancestral evolutionary branch of the genus 
Odorrana. Therefore, we believe that the ancestral distribution 
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of Odorrana may be the region south of Guizhou and northwest 
of Guangzi. However, this preliminary speculation needs to be 
supported by data from additional species and nuclear genes. 

Biodiversity conservation in southwestern China is a priority 
of the Chinese government (Ministry of Environmental Protection, 
2015). Biodiversity conservation programs in this region play an 
important role in maintaining the stability of mountain ecosystems 
as well as protecting biodiversity. In the past three years alone, 
16 new amphibian species have been described from Guizhou 
Province, China (Zhang et al, 2017; Li et al, 2018a, b; Li et al, 2019a, 
b; Li et al, 2020a, b; Lyu et al, 2019; Wang et al, 2019; Luo et al, 
2020; Liu et al, 2020; Lyu et al, 2020; Su et al, 2020; Wei et al, 2020; 
Wang et al, 2020; Cheng et al, 2021). The discovery of these new 
species suggests that amphibian species diversity in this region 
is severely underestimated. In the context of global warming, 
there is an urgent need for a comprehensive, systematic, and 
in-depth survey of the impacts of climate change on terrestrial 
vertebrates to provide a basis for scientific decisions regarding 
amphibian conservation (IPCC, 2014). 
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Appendix 


Specimens examined 

Odorrana lipuensis (n = 8): China: Guangxi: Lipu County (type locality): four males: NHMG 1306001, NHMG 1306002, NHMG 
20140702, and NHMG 20140703; four females: NHMG 1303019, NHMG 1303018, NHMG 1306003, and NHMG 20140701. 
Odorrana kweichowensis (n = 35): China: Guizhou Province: Jinsha County (type locality). 18 males: GZNU 20170717008-010, 
GZNU 20170717014, GZNU 20170717016-019, GZNU 20170717021-023, GZNU 20170717025, GZNU 20170717026, GZNU 
20170717028-7030, GZNU 20170718003, and GZNU 2017071800; 17 females: GZNU 20170717001-007, GZNU 20170718001-002, 
GZNU 20170725001-007, and GZNU 20170725013. 

Odorrana wuchuanensis (n = 17): China: Guizhou Province: Libo County (topotype locality): 12 males: GZNU 20160729001-006, 
GZNU 20160729008, GZNU 20160806001-005; five females: GZNU 20160729007, GZNU 20160809001, GZNU 20160809003, 
GZNU 20160809006, and GZNU 20160809011. 


